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Abstract We study the existence of singular solutions to the equation —div{\Du\^ ^Du) = |it|'' 
under the form u{r, 9) = r-''^uj{e), r > 0,6 e S'^'K We prove the existence of an exponent q below 
which no positive solutions can exist. If the dimension is 2 we use a dynamical system approach to 
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1 Introduction 

The study of isoslated singularities of solutions of quasilinear equations started with the 
celebrated works of Serrin [20] [21] dealing with expressions such as 

divA{{x,u,Du)) + B{x,u,Du) ^0 (1.1) 

where A and B are respectively vector valued and real valued Caratheodory functions sat- 
isfying the same power p-growth with p > 1. One of the main results of these works stated 
that the type of singularities is dictated by the diffusion operator A. Later on the particular 
cases of superlinear semilinear elliptic equations was considered, either with an absorption 

-Au+lul"^'^ u^O (1.2) 

[5] , [24] , or with a source reaction 

Au + u« = (1.3) 

[17], [10], [2], and in all cases q > I. One of the main facts of these studies relied in the 
existence of critical thresholds where the interaction of the diffusion and the reaction terms 
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could create unexpected phenomena. As a natural generalisation, the same analysis was 
carried on for 

-div(|£>M|^"^£>u) + |m|«"^m = (1.4) 

[9], and 

div {\Duf~'^ Dv)j + u« = (1.5) 

[22] , in the range Q < p — 1 < q. In all these works, the radial explicit solutions, whenever 
they exist, played a key role. 

Similarly, the study of the boundary behaviour of solutions of quasilinear equations, has 
a natural starting point in the description of their isolated singularities on the boundary. 
Besides the historical results of Fatou, Herglotz and Doob on the boundary trace of posi 
five harmonic and super harmonic functions, equations of types (1.2 ), (1.3 ) and (1.4 ) have 
alredy been considered ([4], [11], [6]). In the present article we consider equations of type (1.5 
). The problem can be stated under the following form: Assmnc is an open subset of 
M^, a € 90 and u e C(f2 \ {a}) fl C^(f2) is a solution of one of the above equations which 
vanishes on 9fi\ {a}, what is the behaviour of u{x) when x — » a. The simplest configuration 
corresponds to f7 = K;'^, and a = (or more gcncraly, if VI is a cone and the singular point 
a its vertex 0) . For such geometry, the key-stone element for describing the behaviour of u 
near is played by separable solutions, whenever they exist. These solutions, which have 
the form 

u{t) = M(r, a) = r-f^uj{(7) r > 0, cr G (1.6) 

have already proved their importance for (1.2 ), (1.3 ) and (1.4 ). It is expected that such 
will be the case for (1.5 ), even if the full theory will be much more difficult to develop 
because of the absence of comparison principle and a priori estimates near a; = 0. It is 
straightforward that, if u is a separable solution of (1.5 ) in R^, 

P = — T ■■= f^i' (1-7) 

q+l-p ^ 

which is positive since q > p — I. Furthermore w is a solution of 

-v. (^[py + 1 vwi')"^'"' - icof-' uj = Ag,p [py + \vujfy^^~' w, (i.s) 

in S^~^, where V' is the covariant gradient on S^~^, V. the divergence operator acting on 
vector fields on S^~^ and 

Xg,p = PgiqPg - N). 
When p = 2, pg = 2/{q - 1) and (1.8 ) becomes 

-A'w- |a;|«"^w = A,,2W, (1.9) 
where A' is the Laplace-Beltrami operator on S^^^ and 

2 / 2q 
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If 5 is a subdomain of ^, equation (1.9 ), considered in S, is the Euler-Lagrange variation 
of the functional 

m = I [I m' + ^V'^ - ^ li^r') da. (1.10) 

For any 1 < q < {N + 1)/{N - 3) (any g > 1 if iV = 2 or 3) this functional satisfies the 
Palais-Smale condition. Furthermore, if Ag^2 < is the first eigenvalue of —A' in 

Wq'^{S)), Ambrosetti-Rabinowitz theorem [1] or Pohozaev fibration method [18], [19] apply 
and yield to the existence of non-trivial positive solutions to (1.9 ) in S vanishing on dS; 
while if Ag,2 > As,2 no such solution exists. 

When p ^ 2, equation (1.8 ) cannot be associated to any functional defined on 
except ii q = Qc = {N{p — 1) + p)/{N — p) (the critical Sobolcv exponent for W^'^, when 
N > p); therefore, finding functions satisfying it is not straightforward. Besides the constant 
solutions which exist as soon as qf3q < N, it is not easy to prove the existence of non-constant 
solutions. As in the case p = 2, it is remarkable to see that existence, or non existence, of 
solutions of (1.8 ) is associated to some spectral problem, although this problem is not 
standard at all: if one looks for the existence of a positive p-harmonic function v in the cone 
C's ~ {(''j"') '■ r > 0, a € S} vanishing on dS, under the form v{r,a) = r^^(f){a), one finds 
that is a positive solution of the so-called spherical p-harmonic spectral equation on S, 
namely 

-v. (^(/32</.2 + I Vc/.]')"^'"' = A (/?2</,2 + \Vct>fy^'~' <^ in S ^ 

(l) = in dS, 

and X = P (Pip — I) + p — N). The difficulty of this problem is two-fold since (3 is unknown 
and (1.11 ) is not the Euler-Lagrange equation of any functional. However, given a smooth 
subdomain S C S^^^, it is proved in [25], following a shooting method duo to Tolksdorff 
[23], that there exists a couple (/3, 4>) = (/?s, 4>s)j where /3s > is unique and (j)s is defined 
up to an homothethy, such that (1.11 ) holds. Denoting 

\s = f3s{l3sip-l)+p-N), 

the couple {(ps, As) is the natural generalization of the first eigenfunction and eigenvalue of 
the Laplace-Beltrami operator in Wq''^{S) since As = As,2 when p = 2. Our first theorem is 
a non-existence which extends the one already mentioned in the case p = 2. 

Theorem 1. Let S C S^~^ be a smooth subdomain. If Pq > Ps there exists no positive 
solution of (1.8 ) in S which vanishes on dS. 

Apart the case p = 2, the existence counterpart of this theorem is not known in arbitrary 
dimension, except \i q = qc in which case (1.5 ) is the Euler-Lagrange equation of the 
functional 

m = I {pij' + ivvi^^' - ^ i^r=+^) da, (1.12) 
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and applications of the already mentioned variational methods lead to an existence result. 



However, when N = 2 the problem of finding solutions of (1.5 ) under the form (1.6 ) 
can be completely solved using dynamical systems methods. In order to point out a richer 
class of phenomena, we shall imbed this problem into a more general class of quasilinear 
equations with a potential, authorizing even the value p = 1. This equation is the following, 



div ( \Duf-^ Du) + lul"-^ u--^ \uf-^ u = 
\ I \x\ 



(1.13) 



in \ {0}, with g > p — 1 > and c e R. If w is a solution under the form (1.6 ), /? is be 



equal to /3q, while w is any 27r-periodic solution of 

2\ P-2)/2 



d 

da 



duj 
da 



duj 
da 



duj 
da 



(p-2)/2 



|a;|«-iw- 



c\uj\P-^LU = 0, 



(1.14) 



where 
If we set 



_i(p-2)g + 2(p-l) 



(1.15) 



(1.16) 



{q+i-pr ' 

then, if c < c^, the only constant solution is the zero function, while if c > c^, there exist 
two other constant solutions ±(c — Cq)^/^'"'"^"^''. Let us denote by f + the set of positive 
solutions of (1.14 ) on S^, £ the set of sign changing solutions and JF = ±5+ U £ the set of 
all nonzero solutions. Our main result which gives the struture of the sets £ and f + is the 
following: 



Theorem 2. Assume p > 1, q > p — 1. Then 
fi) 



£= [J {ivki. + ij) : ij e S'} 



(1.17) 



fceN, 



in which expression ojk is a function with least period 2TT/k, and kg = 1 if c > Cg, or kg is 
the smallest positive integer such that kg > Mg, where 



Mg = 



Jo I 



l + {p-l) tan^ e 



(1.18) 



/?^(p - 1) tan^ e + Cg- ccosP-2 6 



d0 



if C< Cg. 

(ii) If c< Cg, f + is empty. IfO<c—Cg < /p, is reduced to the constant function 



- r-_^V(9+l-p) 



If c — Cg > ^/p, contains the constant function (c — c, 



- rg)l/(</+l-P) 
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and the set 

fe+ 



U K(. + V'):V'e5^}, (1.19) 



ken, 
fe=i 



where uj^ is a non-constant positive function with least period 2'K/k, and k'^ is the largest 
integer smaller than {pPg~^{c — Cg))^/^ 



Since separable solutions of (1.5 ) defined in a cone Cs and vanishing on dCs are asso- 
ciated to elements of f , we can prove the existence counterpart of Theorem 1 in dimension 2. 

Corollary 1. Let N = 2 and S be any angular sector of . Then there exists a positive 
solution of (1.8 ) vanishing at the two end points of S if and only if Pq < Ps- Further- 
more this solution is unique. In particular, existence holds for any sector if p < 2 and 
q>2{p-l)/{2-p). 

The case p = 1 appears as a limiting case of the preceding one. In that case we observe 
that li is a positive solution of (1.13 ) if and only if = is a solution of the same equation 
relative to g = 1, 

Avv(\Dv\~'^ Dv\+v- = Q. (1.20) 

The initial case c = is easily treated, but the case c ^ Q, that we shall analyse in full 
generality, is much richer and delicate and shows a large variety of solutions depending on 
various parameters. 

Theorem 3. Assume p=l and q > 0. Then 

(i) If c ^ 0, or c = and q> 1, £ is empty. If c = and q <1, £ = {a;o(. + il>) : ip G 5^} , 

where a a'o(cr) := 2^^'^ Isino-I*^^"^-*^^ sin cr is a solution of (1.14 )■ 

(a) If c < —1, is empty. If — 1 < c < Q, £^ is reduced to the constant function (c+l)^/'. 

//c>0, 

f+ = {(c+i)i/''}u U K(. + ^):V'e^^}, 



k—k\ 

in which expression uj^ is a positive function with least period 2'K/k, k2 is the largest integer 

strictly smaller than (c+1)^^^ and ki is the smallest integer greater than — y cose+2c '^^- 
Finally, if c = 0, 

, f if 9 > 1 

£+ = {l}U U {a;+(. + V):Ve5i}u ^ 

Ke(o,i) I I'^o 



i(0,l) 

where the functions and Wq are explicitely given by 

1/9 



iv^= (^Vl-K^ sin^ (7 - if cos a) ^ and = (2 |sincr|)^/« Va G 5^ 
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A striking phenomenon is the existence of a 2-paranieter family of solutions when c = 0. 



Our paper is organized as follows: 1- Introduction. 2- The N-dimensional case. 3- The 2-dim 
dynamical system. 4- The case p> 1. 5- The case p= 1. 

2 The N-dimensional case 

2.1 The spherical p- harmonic spectral problem 

If p > 1, (3 > and A e IR we denote by 1/3, x the operator defined on C^{S'^~^) by 

cp ^ T/3,aM = -v. ({/3V' + |VV|')^^-'^/'VV) - A(/3V' + |VV|')(^-')/V- (2.1) 

Let q>p— 1>0, S he a, smooth connected domain on and Cg the cone with vertex 

generated by S. If u is a positive solutions of 

-diy{\Du\P-'^Du) = u", (2.2) 

in Cs \ {(0)} vanishing on dCs \ {(0)}, under the form 

u{r,a) =r-^Lo{a), (2.3) 

then /3 = p/{q + I — p) := fiq and w solves 

a; = on dS, 

where 

A,,j,=/3,(9/3,-7V). 

We denote by (3s the exponent corresponding to the first spherical singular p-harmonic 
function and by 0s the corresponding function. Thus f3s > and u{r,a) = r^'^^' (ps{o') is 
p-harmonic in Cs \ {(0)} and vanishes on dCs \ {(0)}. Furthermore ^ = ^s > and satisfies 

f 2:fe,AsM =0 in s 
[ (f) = on dS, 

where 

Xs=0sWs{p-l)+p-N). 

We recall that {Ps, <Ps) is unique up to an homothety upon ^. Furthermore ^s is positive 
in S, d(j)s/dv <Q on dS and 

S'^S,S'i^S^ Ps' > (3s. 
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2.2 Non-existence 

Proof of Theorem 1. We put 

^ = ^ and r? = <^%. 

Then ^ > 1 and 

+ {p- i){e - i)ef-vr'^^'"'^"'(/3|'^s + 1 v'<^5|')(^-2)/' I V'-^sl' 

Using (2.5 ) with = ^5, we derive 

=-(p-l)^^-'(^-l)4^"'^^'"'^"'(/3|0| + |V>s|'r/' in 5. (2.6) 

Because a; is a nonnegative nontrivial solution of (2.4 ), it is nonpositive in S. Furthermore 
doj/dv < on dS. Therefore we can choose (j)s as the maximal positive solution of (2.5 ) 
such that T] <w. If 6 > 1 there exists a* G S such that 

w(f7*) = rj{a*) > and uj{a) > r]{a) Va e S. (2.7) 

If 9 — 1, the graphs of co and r] could be tangent only on dS. This means that either (2.7 ) 
holds, or there exists a G dS such that 

du){a)/du = dr]ia)/ dp <0 and Lo{a) <Tj{(7) \/a G S. (2.8) 

Let tp = LU — ri and we first consider the case where (2.7 ) holds. Let g = (gij) be the metric 
tensor on S^~^. We recall the following expressions in local coordinates aj around a*, 

for any (f G C^{S), and 



V'.X = 




for any vector field X G C'^{TS^-'^), if we lower the indices by setting = '^g^'Xi. We 

i 

derive from the mean value theorem 

iPy + |V'c.|^)(-^)/^^ - + |V',y|^)(-^)/^g- = E4^^ + b^co - v), 
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where 



, (p-4)/2 



b' = {p- 2) + t{u; - v)f + + t{co - v))f) 

d{r] + t{uj-r])) 



x{r] + t{uj - m))- 



dcFi 



and 



(p - 2) (plii^ + t(uj - ij)f + |V'(r, + tiu ~ r;))!' 



(p-4)/2 



dui ^ dak 

k 

. / „\ (p-2)/2 

+51 [filir^ + tiu; - rjjf + \V{v + t{co - v))f) 
Since the graph of r] and u are tangent at a* , 

r]{a*) = oj{a*) = Pq > and Vrjia*) = Vuj{a*) = Q. 

Thus 



and 



a]{a*) = (filP^ + 101')^'"'^^' ( 5l{(3lP^ + \Q\^) + (p - 2)g,^5^'=0fe ) . 



Now 



, ,i2xp_i5a; 2^2 , |w„|2.p_i a?] 



\9\ 



dai 



da, 



— T]) 
dai 



where the Cj are continuous functions and 
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The matrix (a* ((7o)) is symmetric, definite and positive since it is the Hessian of the strictly 
convex function 

p/2 

x = {Xi,.. .,Xn-i) ^ - [pS + \xfy^^ = - ( ^0 + E a'^^j^k 

Therefore (aj) has the same property in some neighborhood of a*, and the same holds true 
with (flj). Finally the function tp = oj — t] is nonnegative, vanishes at a* and satifies 



dtp 



dtp 



ddi 



Then i/» = in a neighborhood of 5. Since 5 is connected, tp is identically which a 
contradiction. 

If (2.8 ) holds, then 9 = 1 and the graphs of rj and co are tangent at a. Proceeding 
as above and using the fact that dij/dv exists and never vanishes on the boundary, we see 
that tp = T] — U! satisfies (2.9 ) with a strongly elliptic operator in a neighborhood Af of a. 
Moreover tp>0mj\f, tp{a) = and dtpldv{a-) = 0. This is a contradiction, which ends the 
proof. □ 

Remark. If p = 2, the proof of non-existence is straightforward by multiplying the equation 
in w by the first eigenfunction (ps and get 

/ {{\s-\,2)^-^'')(psda = Q, 
Js 

a contradiction since Ag < Ag,2- 
2.3 Existence results 

Let us consider the case q = Qc = {N'{p — 1) + p)/{N — p) {N > p > 1), and let S be any 
smooth subdomain of S^~^. Since in that case Ag,p = — /3q^, the research of solutions of (1.5 
) under the form (1.6 ) vanishing on dCs leads to 

f '^0.c,-0L M - I'^l'^"' Lo = inS 
\ LU = in dS, 

where (3q^ = N/p — 1. This equation is the Euler-Lagrange variation of the functional J 
defined on Wo'^(5) by 

J{tP) = g (Pltp' + \Vtpff" - ^ |^r=+^) da. (2.11) 

Theorem 2.1 Problem (2.10 ) admits a positive solution. 
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Proof. Clearly the functional is well defined on Wq'^{S) since Qc is smaller than the Sobolev 
exponent p*j^_i for W^'^ in dimension N-1. For any G Wq''^{S), liiat^ao J (tip) = — oo. 
Furthermore there exist S > and e > such that J{ip) > e for any ip e Wq''^{S) such that 
ll^llwi.p = ^- Assume now that {ipn} is a sequence of Wq'^{S) such that J^ipn) — * o. and 
||-DJ(f/;„)||^_i,p' ^ as n ^ oo. Then 

a it follows 

(/3,lV'n + I VVnl')"^' da ^ p(5e + l)a/{qc + 1 - p). 

Therefore {V'n} remains bounded in L''^+^(S'), and relatively compact in L'^{S), for any 1 < 

r < Qc + l. Multiplying the equation _DJ(?A„) - e„ by Tfe^eCV-'n) where 6* € (1, {p*j^_i - l)/<?c), 
fc > and Tk^0{r) = sgnmin{|r| , k} and using standard bootstrap arguments yields to the 
boundedness of {ipn} in L°°{S). Combining this fact with the compactness of in L'^{S), we 
derive the compactness in any any L^, for s < oo. Therefore {ipn} is relatively compact in 
Wo'P{S). This means that J satisfies the Palais-Smale condition. □ 



Then 

Since J(ipn) 



3 The 2-dim dynamical system 
3.1 Extension of the data 

Due to possible applications and similarly to what is done in the semilinear case p = 2 
(see [3], [7], [8]), we shall consider the existence problem for 27r-periodic solutions of a more 
general quasilinear equation than (1.14 ), 



d_ 
da 



■2, ,2 



duj\ 



2\ 



duj 
da 



dio 
da 



p/2-l 



+ g{w) - c\w\P~'^uj = 0, 



where X,l3,c are real parameters, with f3 > 0, and g G 
satisfies 



n 



(3.1) 

{0}) is odd and 



lim o(s)/s'' = 1, 

«— >o+ 



lim g[s)/s^ ^ = oo 



d_ 
ds 



(5(s)/|sr') >0on (0,oo), (3.2) 



with g > p — 1 > 0. In fact we can easily reduce the problem to a simpler form, and 
particularly in the case p = I, where the equation has a remarkable homogeneity property. 
The next statement is a straightforward computation which transforms the equation satisfied 
by u) into two more canonic forms. 
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Lemma 3.1 Let oj be a solution of (3.1 ). 
(i) Assume p > 1. If we set 

T = p(7, a;((T) =/?f/(«+i-P'u;(r) andw' = ^, (3.3) 

dr 

then w satisfies 

^ ^ ^/2)P/2-l _ (^2 ^ «;'2)f/2-l ^ ^ ^(^) _ rf|^|p-2^ ^ (3 4) 

where 

b=^^ d=-^^ /(s)=r^''/(«+i-P)5(/J^'/(''+i-f)s). (3.5) 

In particular f satisfies the same assumptions (3.2 ) as g. 

(a) Assume p > 1. If on any open interval I C (0, 27r) where u}{a) ^ 0, we set 

T = f3qa, anda;(a) = (/3g)i/« |w(T)|^/«~^w;(r), (3.6) 

then w satisfies ( 3.4 ) on I, with 

b=-X/p\ d = c/pq, /i(s) = /J-''5((/3gs)V9). (3.7) 

Furthermore fi satisfies the assumptions (3.2 ) with q= 1, i.e. 

lim fi(s)/s = 1, lim /i(s) = oo, fUs) > on (O.oo) . (3.8) 

Due to this result, the changes of variables (3.3 ) and (3.6 ) reduce the problem to the 
study both of existence of periodic solutions of equation (3.4 ), and to characterizing the 
period function of these solutions, in the range q> p — 1 if p > 1, and g > if p = 1. 

3.2 Reduction to dynamical systems 

We re-write (3.4 ) as the system, 
w' = F{w, y)=y 

, _ bw^ + {b + 2- p)w y'^ - (/(w) - d\w\P-^w){w^ + y'^f-p/^ (3.9) 

+ {p- l)j/2 

and we denote by h the odd function defined on R by 

h(s) = J ^ I I ^ (3.10) 

ifs = 0. 
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l{ b + d < 0, (3.9 ) has no non-trivial stationary point, while if 6 + d > 0, it admits the two 
stationary points ±Po, with Pq — (a, 0) and a = h~^{b + d). Furthermore Pq is a center since 
the linearized system at Pq is given by the matrix 

1 

-ah'{a) 

System (3.9 ) is clearly singular at (0,0). Furthermore it could singular be along the line 
w = Q \i p = 1, \i q < 1, and lip <2 and d^Q. Actually, for p > 1 it is not singular at any 
points (0, a) with a j^O. This can be checked as follows: consider the Cauchy problem 

w" = G{w,w'), te{-d,5) 

wiO)=0,w'{0) = a, ^'^ > 

and let w be any local solution; since near (0, a), G is continuous with respect to w and 
with respect to y, w is C^; because a ^0,t can be expressed locally in terms of w. Defining 
w'{t) = p{w), then p is near 0, p{0) = 1 and satisfies 

dp G{w,p) 
dw p ' 

with J{w,p) — G{w,p)/p. Clearly is with respect to p and continuous with respect to w, 
thus one gets local uniqueness of p. and then the local uniqueness of problem w'{t) = p{'w{t)), 
w{0) = 1, since p is of class C^. 

The phase plane of the system (3.9 ) is equivariant under symmetries with respect to 
the two axes of coordinates, because F is even with respect to w and odd with respect to y, 
and G is odd with respect to w and even with respect to y. Thus from now we can restrict 
the study to the first quadrant 

Q\{(0,0)}, where Q = (0, oo) X (0, oo) , 

where, in particular, w > 0. Due to the symmetries, in the case p > 1,, any trajectory which 
meets the two axes in finite times r, r + T is a closed orbit of period 4T. 

Remark. It is useful to introduce the slope ^ — w' /w, (or a function of the slope) as a new 
variable. This was first used for p > 1 in [16] for the homogeneous problem 

w'Y^-' w') -b{w-'+ w'Y^-' w = 0. 
In that case the function ^ satisfies 

A ((1 + ^Y'-' ^) = -iiP- i)e -b){i+ ef'-' , 

for w > 0, and this equation is completely integrable in terms of u = (l + ^^Y^^ ^ ^■ 



12 



By using polar coordinates in Q 

{w, y) = [p cos 6,psm9), p>0, 6 € (0, 7r/2), 
we transform (3.9 ) into 

b-{p-l) tan2 e+{d- h{p cos 9)) cos^'-^ 9 

l + (p-l)tan20 (3.12) 
p' = p{l+9')i&n9. 

Equivalently, if we introduce the slope ^ = tan^ e (0, oo), and set 

u = (f){Cj = cos^-P9sm9, (t){0 = (1 + ^^)^P"^^/^^, (3.13) 

then (^'(^) = (1 + ^2-)(p-4)/2^-[^ + (p — 1)^^); thus (j) is strictly increasing: from (0,00) into 
(0, 00) when p> 1, and from (0, 00) into (0, 1) when p = l. Defining 

ip = r\ and E{0 = {{p-l)e-b){l + er^^-\ (3.14) 

we obtain 

{w' = wip(u), 
(3.15) 
u' = -E{ip{u)) - h{w) + d. 

This system is still singular on the line w = if h ^ {[0, 00)) near 0. In the sequel we set 

u 

*(u) = j ip{s)ds. (3.16) 


Noticing that 

E'iO = {piP - + 2(p - 1) - (P - 2)6) (1 + ^2)(P-2)/2^^ (317) 

we derive that E is increasing on (0, 00) when {p — 2)b < 2{p— 1). When (p — 2)6 > 2(p— 1), 
E is decreasing on (0, rj) and then increasing, where r] is defined by 

p{p-l)r]^ = (p-2)b-2{p-l), (3.18) 

and 

minE = EM = (iP-m^P-DV^' 
In the case of initial problem (1.14 ), E is increasing. 

Remark. If p > 1, system (3.9 ) is singular at (0, 0). If we replace the assumption lims_>o+ f{s) /s^ = 
1, by the stronger one 

hm f'{s)/s'i-' = q, (3.20) 

s—>0+ 
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we can transform system (3.15 ) in (0, oo) x M in a system of the same type, but without 
singularity: this is obtained by performing the substitution v = w'^~^^~p. Then 

/ v' = {q + 1 - p)vip{u) 

\ u' = -E{ip{u))-hiv)+d, ^ ' 

where v i— » h{v) = e C^([0, 1) . In particular, if f(w) = \w\'^~^w, we find 

( v' = {q + 1 - p)v(p{u) , . 

\ u' = -E{ip{u))~v + d. ^'^ > 

Remark. In the case f{w) = \w\''~^w, we can differentiate the equation relative to u' and 
obtain that u satisfies the following equation 

u" = B{ip{u))u' + {q + l- p){E{(p{u)) - d)(p{u), (3.23) 

where E is given above, and 

o..^ (p-2)6 + g-3(p-l) + (g + l-2p)(p-l)^% 

m = i + {p~i)e ^- ^ ^ ^ 

Notice that equation (3.23 ) has no singularity for p > 1. 

4 The case p > 1 

4.1 Existence of a first integral 

A natural question is to sec if equation (3.4 ) admits a variational structure. When p = 2, 
it is the case, for any b and d. Since (3.4 ) takes the form 

w" -{b + d)w + f{w) = 0, 

it is the Euler equation of the functional 

n2{w,w') = ^ + {b+d)^-J^{w) 

pw 

where ^(w) = / f{s)ds. Thus the function w''^ = {b + d)w'^ — 2J^{w) is constant along the 
Jo 

trajectory. When p 7^ 2,p > 1, we find that a first integral exists only in the case 6 = 1. In 
such a case (3.4 ) is the Euler equation of the functional 

H{w, w') = ^ '- + rfJ— L - T{w). 

p p 

Therefore, the associated Painleve integral 

Viw, w') ^-(w^+ w'Y^^' i(P - - «^') - + (4-1) 

p P 
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is constant along the trajectories. Using the function E introduced at (3.14 ), then (4.1 ) is 
equivalent to 

E(^)=i5M„» = <,-pi:±|M (4.2) 

for w > 0. Hence E is increasing on (0,oo) from — fe = —1 to +oo. 

In the general case, we cannot use a first integral for studying the periodicity properties 
of the solutions, while it was the main tool in [3] for p = 2. This is the reason for which 
wc arc lead to use phase plane techniques. Notice that, for the initial problem (1.14 ), the 
value 6=1 corresponds to the case p <2 and q = {3p — 2)/(2 — p)). 

4.2 Description of the solutions 

In this section we describe in full details the trajectories of system (3.9 ) in the phase plane 
{w,y). Notice that the system can be singular on the axis w = 0. 

Proposition 4.1 Assume p > 1. Then all the orbits of system (3.9 ) are hounded. Any 

trajectory T[p] issued from, a point P in Q is 

(i) either a closed orbit surrounding (0,0), 

(ii) or, if b + d > 0, a closed orbit surrounding Pq but not (0, 0), 

(Hi) or an homoclinic orbit defined on M, starting from (0,0) with initial slope 

lim — — = m 

t^-oo w{t) 

where m is defined E{m) = d, and ending at (0, 0) with 

lim — — = —m. 

t->cx) w{t) 

Proof. We recall that E and u are defined by (3.13 ) and (3.14 ), by using polar coordinates 
{p, 0) in the {w, T/)-planc. 

First look at the vector field on the boundary of Q. At any point (0, a) with cr > 0, it is 
given by {a, 0), thus it is transverse and inward. At any point {w, 0) with w > 0, it is given 
by (0, w{b + d — h{w)) . Thus it is transverse and outward whenever b + d<OoTb + d>0 
and iu > a, and inward whenever b + d> and w < a. 

Consider any solution {w,y) of the system, such that P = (w(0),y(0)) G Q, and let 
(''"ii T2) be its maximal interval existence in Q. At any point t where u'(t) = and u(r) > 0, 
there holds u"{t) = —h'{w)w(p{u) < from (3.15 ). Thus if r exists, it is unique, and it is 
a maximum for u. 

Since w' = y > 0, w has the limits £2 € (0, 00] as t t r2 and £1 S [0,oo) as r | ti. 
Therefore u is strictly monotonous near n, and T2 thus it has limits Wi, U2 € [0, 00] , in other 
words 6 has limits 61, 02 € [0, 7r/2] 

(i) Let us go forward in time. On any interval where u is increasing, one has E{(p{u)) < d, 
thus u is bounded and, consequently, U2 is finite. If ^2 = oOj then 0'{t) — > — 00, as r t T2; 
by (3.13 ), p is decreasing, thus it is bounded, which is contradictory; thus £2 is finite. If 
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U2 > then {£2,^2<f (^2)) is stationary, which is impossible. Thus u is decreasing to 0, and 
the trajectory converges to (£2, 0). If 6 + > and £2 = a, u' tends to from (3.15 ), and 

u" = —{E o ip)' (u)u — h' (w)w(p{u) = —h'{a)a(f{u){l + o(l)); 

therefore u" < near T2, which is impossible. Finally, either 6 + d<0, or6 + d>0 and 
id > a, and T2 is finite, the trajectory leaves Q transversally at T2. 

(ii) Next let us go backward in time. 

• Suppose ui = 0. Clearly the trajectory converges to {£1, 0); then necessarily b + d > 
and £1 < a, thus < a as above. The trajectory enters Q transversally at T2, and from the 
symmetries it is a closed orbit surrounding only the stationnary point Pq. 

• Next, suppose ui = 00. It means that 9 tends to 7r/2. Then from (3.12 ), 6' tends to 
1, thus Ti is finite, 

7r/2 - = (r - ri)(l + o(l)) , tan0 = (r - ti)-^{1 + o(l)), 

and ^ 

(p - 1)(t - ri)-i^ = {b+l + {d- h{p cos 6)) cosP-^e)(l + o(l). 
P 

If p > 2, then p' /p = 0((r — Ti)); if p < 2 then p'/p = 0((t — ri)^"^). In any case, \np case 
is intcgrablc, thus p has a finite limit ^ > 0. Then the trajectory enters Q transversally at 
Ti and from the symmcrics it is a closed orbit surrounding (0,0). From the considerations 
in § 3-2, for any y > there exist such an orbit, and it is unique. Moreover in Q the slope 
w' /w = ^ = ip{u) is decreasing from 00 to 0; indeed it decreases near n and T2 and can only 
have a maximal point. 

• At end, suppose < ui < 00. If fi > 0, then {£i,£iLp (ui)) is stationary, which is 
impossible. Thus {y,w) converges to (0,0). And w'/w tends to f{ui) , thus n = —00. And 
u' converges to d ^ E{ip[ui)), thus tanfl — (p(u) has a limit m > such that E{m) = d. 
From the symmetries the trajectory is homoclinic and the solution w is defined on R. □ 

The next theorem studies the precise behaviour of solutions according to the sign ofh+d. 

Theorem 4.2 Assume p> 1 and consider system (3.9 ) in the {w,y)-plane. 

(i) Assume b + d > 0. Then there exists a unique homoclinic trajectory Ti. starting from (0, 0) 
in Q with initial slope mj, = E~^{d) ( TOq = \Jb/{p — 1) if d = 0), ending at (0,0) with 
the slope —m,d, and surrounding Pq- Up to the stationary points, the other orbits are closed, 
and either they surround only one of the points Pq or —Pq, in the domain delimitated by 
Ti., corresponding to solutions w of constant sign, or they are exterior to zLH and surround 
(0,0) and ±Po, corresponding to sign changing solutions w. 

(ii) Assume b + d < 0. Then 

• if {P ~ 2)^ < 2(p — 1), or [{p — 2)6 > 2(p — 1) and d < -E'(ry)] , there is no homoclinic 
trajectory. 

• */ [ip ^ 2)6 > 2{p — 1) and E{ri) < d < —b] ; then denoting by mi.^ < m,2.d the two 
positive roots of equation E[m) = d, there exist infinitely many homoclinic trajectories Hi 
starting from (0,0 in Q with the initial slope mi,d and ending at (0,0) with the final slope 
— mi_(i, and a unique homoclinic trajectory TI2 starting from (0,0) in Q with initial slope 
777-2, d and ending at (0, 0) with final slope —m2,d- 
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Proof, (i) Case b + d > 0. Then the equation E{m) = d has a unique positive solution 
m = E~^{d); and w' /w tends to m; thus the trajectory starts from (0,0) with a slope m. 
Then for any P G Q, the trajectory 7[p] passing through P meets the axis j/ = after P at 
some point (/x, 0) with fj, > a. Denote 

U={PgQ: T[p] n {(0, ct) : CT > 0} ^ 0} , V= {P e Q : T^P] n {{n, 0) : < m < a} 7^ 0} • 

(4.3) 

Then either P gU and the trajectory is a closed orbit surrounding (0, 0) and ±Po, and in Q. 
Or P G V and the trajectory is a closed orbit surrounding only Pq. Or 7[p] is an homoclinic 
orbit Ti. starting from (0,0) with the slope m, where m is the unique solution of equation 
E{m) = d (such that m > if £^ is not monotone, see (3.17 )). Next U and V are open, 
since the vector field is transverse on the axes, thus UUV ^ Q. This shows the existence of 
such an orbit Ti. 

(ii) Case h + d<Q. 

• Either b + d < and E is increasing, or E has a minimum at 77 and d < E{r]). In 
such a case equation E{m) = d has no solution, and there is no homoclinic orbit. Or E is 
increasing and 6 + d = 0; then E{(p{u)) > —b = d, thus u' < 0, thus u cannot tend to 0, and 
the same conclusion holds. 

• Or E has a minimum at rj and E{r]) < d < —b. In that case the equation E{m) = d has 
two roots mi, 1712 such that < mi < rj < 1712 < ruh, where m{b) is defined by E{mb) = —b. 
Any trajectory 7[p] such that P E U (sec (4.3 ) for the definition) satisfies u' < 0, it 
means h{w) > d — E{(p(u)) and the range of u is (0,oo), therefore there exists t such that 
(^(u)(r) = T], hence /i(w(r)) > d — E{r]) and y(r) = r]w{r). Next consider any trajectory 
7|pj starting from P = {w, rjw) such that h{w) < d — E{r]). It cannot be a trajectory of the 

preceding type, thus {y,w) (0,0) as r — > n, and 9 tends to 61, with tan^i = mi or 1712; 
moreover u'{0) > 0, and w' < near T2, thus there exists a unique t > such that u'(r) = 0; 
then m' > in (ti, t), therefore tan 9i < rj, and finally tan^i = mi. Consequently there exist 
infinitely many such trajectories Tii, with initial slope mi. Next fix one trajectory ^p^,] such 

that h{wo) < d — E{r]). Let TZ be the subdomain of Q delimitated by ^^Pq] and '^[(0,1)] and 

V={P &n: T^p] n {{w, r]w) : Q< w < Wo} ^ $] . 

The set V is open because the intersection with the line y = rjw for w < wo is transverse 
since at the intersection point, h{'w) < d — E{r]), thus u' > 0, and y/w = ip{u) = t], and 

^ = ^{u) + ^^u' >r,= —. 
y (p{u) w 

Then {U r]TZ)UV ^TZ. Then there exists at least a trajectory Hi^* starting from (0, 0) with 
initial slope m2. 

(iii) Uniqueness of Ti and H2- Let m = mo or m2,d- Suppose that system (3.9 ) has 
two solutions {wi,yi), (w2, 2/2) defined near —00, such that Wi > and Wi(r) tends to and 
yi{T)/wi{T) tend to m as r | —00. Then the system (3.15 ) has two local solutions (^1,^1), 
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{w2,U2) such that ip{ui) tends to m at — oo. Then > locally and one can express Ui as 
a function of Wi. Then at the same point w, 

w — — — = wip{Ui)— = -E{ip{ui)) - h{w) + d, 

^d(*(n2^-«^ = -{E{ipiu2)) - E{v{u,)) = E'{^{u*))^'{u*){u2 - ui) 

for some u* between U\ and 112, and E'{ip{u*)) = E'{m){l + o(l)); and E'{m) > 0. Then for 
small w 

aw 

which implies that {^{u2) — is decreasing, with limit at 0. Therefore ^'(^2) = 

^{ui), thus U2 = Ui near —00; but from (3.15 ), h{wi) = h{'W2), and since h is one to one, it 
foUowswi = W2 near —00. The global uniqueness follows, since the system is regular except 
at (0,0). All the trajectories are described. □ 

Remark. Under the assumption (3.20 ), existence and uniqueness of Ti and W2 can be ob- 
tained in a more direct way whenever d 7^ E{ri). Indeed the system (3.21 ) relative to 
{v,u) is regular, with stationary points (0,0), (0, ±iyj~^(m)), where m = mo, mi or m2 
and also (±a, 0) if & + > 0. The linearized system at (0, ip~^{m)) is given by the matrix 
^m{q + 1 - p) ^0^^^ ^ ^.^^ ^^^^ ^ _ ^^^^2 _ ^2^^^^ ^ _ ^ ^ ^^^^^ 

then it is a source, and we find again the existence of an infinity of solutions. If m = md or 

m — m2.d, then K{m) < 0, thus this point is a saddle point. Then in the phase plane (w, u), 
there exists precisely one trajectory defined near —00, such that v > and converging to 
(0, m) at —00, and u/v converges to 0. □ 

Remark. Suppose /(w) = then we can study the critical case {p — 2)6 > 2{p — 1) 

and E{ri) = d : there exist infinitely many homoclinic trajectories Hi starting from (0, 0) 
in Q with an infinite initial slope and ending at (0, 0) with an infinite slope, and a imiquc 
homoclinic trajectory Ti.2 starting from (0, 0) in Q with the initial slope 77 and ending at 
(0,0) with the slope —77. Indeed using system (3.22 ) and setting u = ip~^{rj) + z, and 
C = (g + 1 — p)r]z + V, it can be written under the form 

Q' = P{C,v), v' = {q + l-p)r,v + Q{Q,v), 

where P and Q both start with quadratic terms. Moreover the quadratic part of P{C,, v) is 
given by P2,oC^ +Pi,iC^ +Po,2W^, where by computation, 

P2,o = -^^^^-^w'^^ir^m+rff^-^y^ < 0. 

q+l-p 

The results follow from the description of sadle-node behaviour given in [13, Theorem 9.1.7]. 

Remark. In the case h=l> —d, we have a representation of the homoclinic trajectory : it 
corresponds to = in (4.2 ). In the case f{w) = \w\''~^w, in terms of u we obtain 

u'= '^'^^^~^ {E{<f{u))-d), 
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which ahows to compute u by a quadrature. 



4.3 Period of the solutions 

First we consider the sign changing solutions 

Theorem 4.3 Assume p> 1. For any v > let 7[(o,j/)] be the trajectory which starts from 
(0, u), and let T{i') he its least period. Then v ^ T(y) is decreasing on (0, oo) . Furthermore 
the range ofT{.) can be computed in the following way. 

(i) Ifb + d< and m E{m) is increasing, or if d < min E, then T(.) decreases from Td 
to 0, where 

T ^iT^' l + (p-l)tan^0 

Jo E{^{u))-d Jo {p-l}ta.nH-b-dcosP-^e ' 

and Td is finite if and only if b + d < 0. If b < = d, then 

{p -1)7 + 1 



^° = '" (7-1)7(7+1) ^ = vW(^). (4.5) 

(a) lfb + d>0 or b + d <0 and d > min£', then T{.) decreases from oo to 0. 

Proof. Step 1. Monotonicity of T. Consider the part of the trajectories '7[(o,i/)] located in 
Q, given by {w^,, y^). We have already shown that u is decreasing with respect to r from oo 
to 0, then E{(p{u)) + /i(wy(u)) — d> and can be expressed in terms of u, and 

/•oo 

^^""^"Vo E{^{u)) + h{w,{u)) - d- ^^-^^ 

Let A > 1. Since the trajectories 7^(o.jy)] and T[(^o,\i^)] have no intersection point, wxi,(u) > 
w^{u) for any u> 0, and h is nondccrcasing, thus T{Xv) < T{v), and T is decreasing. 

Step 2. Behaviour near oo. Let Vn > 1, such that limz/„ = oo. Observe that for fixed m, for 
any integer n > 1, there exists a unique Vn > (depending on u), such that Wv^{u) = n; 
let i^n = max(t'„,n). Then h{w0^{u)) > h{n), thus h{wp^{u)) converges to oo; since v i-^- 
h{wv{u)) is nondecreasing then h{w,^^{u)) converges to oo, and T{vn) converges to 0, using 
the Beppo-Levi theorem. 

Step 3. Behaviour near 0. 

• First assume b + d < 0, and E is increasing, oy d < E{rf). Then all the orbits are of the 
type 7j(o,;/)]. Let Vn G (0, 1) , such that limi/„ = 0. For fixed u and any integer n>\, there 
exists a unique > (depending on u), such that Wj^^{u) = 1/n; let i>n = min(j/„, 1/n). 
Then h{wi,^{u) < h{l/n), thus h{wi,^{u)) converges to 0, and again h(wv„{u)) converges to 
0. Then T{vn) converges to Td given by (4.4 ), using the Beppo-Levi theorem. If 6 + d < 0, 
then Td is finite: indeed near oo, E{(p{u)) — [p — 1)uP/(p~^)(1 + o(l)); if E is increasing, 
then E{ip{u)) - d> -{b + d) > d < E{r]), then E{ip{u)) -d> E{r]) -d>0. 

If 6 + = and E is increasing, then = oo: indeed near 0, 

E{ifi{u)) -d = u2(£;"(o)/2 + o(l)) 
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and 

^"(0) = 2(p - 1) - (p - 2)6 if (p - 2)6 = 2(p - 1). 

Therefore 

E{y,iu)) -d= (pip - 1)/A)u\l + 0(1)). 

In all the cases the integral (4.6 ) giving T is divergent. 
When b < = d, one can compute To : 

To ^ r du ^ rnod^ ^ r i + {p-i)e 
4 io io ^(0 h {\h\ + {p-i)e){i+e) ^ 

= ^ + (^- 2^ r ds _ ^ 1-(P-1)7^ 

2 >-l ^Vo (72 + s2)(l + s2) 2^ ^(p- 1)7(7 + 1)^- 
Hence (4.5 ) holds. 

• Next assume d > E{ri). Considering Un as above, for any fixed u such that ip{u) > TO2, 
there exists a unique > (depending on u), such that Wip„(u) = 1/n. As above, 

du f°° du 



v-^(m2)E{v(u)) +h{wpj - d £((/?(«))- d 

since E' {1112) is finite. As a consequence, T{vn) tends to 00. If = E(ri), the same proof 
still works with m2 replaced by rj : the integral is still divergent because the denominator is 
of order 2 in u — y~^(r/), as, near 0, there holds 

EMuj) - d = \e" (r;) {^{u) ri)\l + o(l)) = \e" (r?) (^(«) - r^f{l + o(l), 

and 

E" (rj) = 2p{p - 1)772(1 + ,,2)(p-2)/2 > o_ 

At last suppose 6 + d > 0; the same proof with 777,2 replaced by m shows that T{i') converges 
to 00 as 1/ tends to 0, since E'{m) at m = E~^{d) is finite. □ 

The monotonicity of the period function is a more general property, since we have the 

following result. 

Proposition 4.4 Let F,G G C^(]R2\(0, 0)) are such that F (resp. G) is odd with respect to 
y (resp. x) and even with respect to x (resp. y), with F{w,y) > in Q. Assume that for 
any {w, y) G Q, and any A > 0, 

d fF(Xw,Xy)\ ^, B (G{\w.\y)\ , , 

' ^ '^^ > > {resp. < 0) and — ^ ' < {resp. > 0). (4.7) 



Assume also that for any a in some interval {cri,a2) (where < ai < 02), the trajectory 
^(0,0-)] of solution of the system 

w' = F{w, y) 



y'^G{w,y) (^-^^ 

passing through (0,cr) (necessarily entering Q since F{0,a) > 0) leaves Q transversally in a 
finite time T{(j)/A at some point (c(a-), 0) (thus G(c(a), 0) < 0). Then (from the symmetries), 
'^(0,0-)] is a closed orbit surrounding (0,0), with period T{a), and a 1— > T{a) is decreasing 
(resp. increasing) on ((7i,(72). 
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Remark. We can notice the condition on F is equivalent to F{Xw, Xy) > XF{w, y) for any 
A > 1. The second condition imphes that for any A > 1, 



G{Xw,Xy) < XG{w,y) {resp.G{Xw, Xy) > XG{w,y)). 



Proof of Proposition 4-4- In polar coordinates (p, 6) in Q, we get 

p' = F cose + G sine, e' = -{G cose - F sine). 

P 

At each point r where ^'(r) = 0, there holds 

, fdG ^ dF . \ . F [dG ^ OF . 

pe"ir) = [^- cose --s.nej , PV) = ^(^cos^--sm< 

But (4.7 ) is equivalent to dF/dp > F/p and dG/dp < G/p (resp >), thus 

P^"{t) < -^—{Gcose - Fsine) = (resp. >) 
pcost' 

In both case 9" has a constant sign. But 9'{0) = -F(0, cr) < and 9'{(7) = G(c(cr),0) < 
thus we get a contradiction by considering the first (resp. the last) point where 9'{t) = 0, 
which satisfies ^"(r) > (resp. < 0). Thus e is decreasing from 7r/2 to 0. Then the curves 
can be represented in function of e by {p {a, e) , e{a)) , and 



a)=4r 
Jo 



1 



H{p{a,e),e) 

with 



H(p, 0) = - iF(p cos e, p sin e) sin 9 — Gip cos 9, p sin e) cos 9) 
P 

Let A > 1. Since the trajectories Tj^o.o-)] T^[(o^\a)]. have no intersection point, then 
p{Xa,e) > p{cr,9) for any £ (0,7r/2); by assumption, for fixed e, the function p 
F{p cos e,p sin 0)/p is nondecreasing (resp. nonincreasing) and p i-^ G{p cos e, p sin 9)/ p is 
decreasing (resp. increasing), thus H{p (Act, 9) ,9) > H{p (ct, 9) , 9), which yields to T(Act) < 
T(ct) (resp. >). This implies that T is decreasing (resp. increasing). □ 

Next we consider the positive solutions w of equation (3.4 ). 

Proposition 4.5 Assume p > 1 and b + d > 0. Consider the trajectories 7[(^ q)]*^ the phase 
plane {w,y) which goes through (/U, 0), for some p, e (0,a) . Let T^{p) he their least period. 
Then 

27r 

lim T+(it) = 00, lim T+(/x) = , 

/'-o t^^a ./ahl{^ 

In particular if f{w) = \w\'^~^w, then lim^^^ T^ip) = 27r/ (g + 1 — p){b + d). 



21 



Proof. We notice that the trajectory ^(^,o)] intersects the Une y = at (/it, 0) and another 
point {g{iJ,),0), with /x < a < g{iJ^), and g is decreasing. 

Stepl. Behaviour near a. When ^, tends to a, then also g{^) tends to a. Indeed for any 
small £ > 0, then g{n) — a < e as soon as /i — a < min(e, a — g~^{a + e)). Since, along such 
a trajectory in Q, ^ = Lp{u) varies from to 0, it has a maximal where u' = 0, thus 
E{^*) = h{w*). When /x tends to a, then h{w*) tends to b, thus ^* tends to E~'^{b) = 0, 
thus also maxj,£7|j^ ^^^j |?y| tends to 0. Using the linearized form of the system at Pq, and 

polar coordinates with center (a,0), w — a + rcosrj, y = \J ah'{a)r sinr], then r tends to 
as /X tends to a, and one finds r/' = ~y/ah'{a) + R/r, where R involves the derivatives of 
G of order 2, which are bounded near the point (a,0), thus i?/r^ is bounded. Therefore rj' 

tends to —yjah'{a), and finally T+(/x) tends to 2-^ / ah' {a) . 

Step2. Behaviour near 0. On the trajectory ^(/j,o)j) the function u is increasing up to a 
maximal value u*{ijl), and then decreasing; moreover u* is a nonincreasing function of fi, 
because two different trajectories have no intersection. Let /x„ G (0, a) , such that lim/i„ = 0. 
For any n there exists /i„ G (0,a) such that the orbit ^(^„,o)]) contains a point above the 
line y = (p~^{m){l — l/n)w, let fin = min(/x„, 1/n). Then u*(/x„) > <p~^{m){l — 1/n), thus 
u*{Hn) tends to m ; then from the Beppo-Levi theorem 

du r°° du 

'(n)E{ipiu)) - d + h{w{iJ,,u)) E{ip{u)) - d + h{w{u)) 

where w is the solution defining H, and this integral is infinite. □ 



f ^ fit J f ^ 

liminfT+(/i) > lim / ^, , ' , , — -, = / 

Ju'(L,)E{ip{u}) - d + h{w[ii,u)) 



Remark. Here the question of the monotonicity of the period is difficult to answer, oven for 
p = 2, where it is solved by using the first integral, see [3]. It is open in the general case. 
More generally, if a dynamical system a center, the description of the period function is still 
a chalenging problem. For example, one can contnict a quadratic dynamical system with a 
center, the associated period function of which is not monotone, and even with at least two 
critical points, see [7] and [8]. 

Remark. In the case 6 = 1, we can compute theoretically the period T+ by using the first 
integral (4.1 ). The stationary point Pq = (/i~^(f ), 0) is obtained for Ka = a^/p — T{a) > 
(in case of a power, Ka ~ (q + 1 — p)/p{q + 1)). The positive solutions correspond to 
trajectories Tr- with K £ (0, Ka) , intersecting the axis y = at points (w;!, 0), {w2, 0) with 
wi < a < W2 defined by w^/p — J^{wi) = K, and the period is given by 

r- = 2- 



rW2 



wE-^{-p^^^y 

Unfortunately, this formula does not allow us to prove the monotonicity of the period func- 
tion for p ^ 2. 

It is remarkable that, in the case f{w) = w, one can solve completely the problem 

in the particular case where 6=1 and q = 2p — 1, using the equation (3.23 ) satisfied by u. 

Proposition 4.6 Suppose that f{w) = \uj\'^~^ w, and 6=1 and q = 2p — 1, p > 1 and 
d + 1 > 0. If p > 2 or d + 1 < 1/ (2 — p), then T"*" is decreasing on (0, a). 
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Proof. Since B{^) = by (3.24 ), equation (3.23 ) turns to 

u" = {q+l- p){E{ip{u) - d)ip{u) = {q+l-p){-{l + d)+ p^.{u))ip{u) 
Henceforth 

^ u"u' = -(1 + d)^'{u)u' + p*(u)*'(u)u', 



q+l-p 

from which expression we derive the first integral, 
1 



q+l-p 



C-U{u)), U = Mo^, M{t) = 2{l + d)t - pt\ (4.9) 



From (4.9 ) the integral curves S in the (u, ■u')-plane are symmetric with respect to the axis 
u' = 0. The times for going from u = to w = w* and from u* to are equal, and u* is 
given by C = M.{^{u*)). The computation of the period is reduced to the part relative to 
the first quadrant. Here we follow the method of [3]: we get 



Jo JUiu*) - UM Jo 



'ds 



^U(u*) - U{r]) ~ Jo ^U{u*)-U{su*)' 
Then 

dT+{u*) ^ {e{u*) - Q{su*))ds 



du* 

and 



dQ{u*) 



= 26' (m*) = U'iu*) - u*U"{u*) 



du* 

In the interval of study, ip{u*) < E~^{d), {Eoip){u) < d from (3.14 ), thus *(«) < {l+d)/p, 
and M is increasing for < f < (1 + d)/p, thus W > 0. Then at any point u, Q'{u) > <S=> 
{U'/uY < 0. Now 

2pu u 
hence {W /u)' = 2X{u)ip(u)ip'{u), with 

X{u) = -{p - 1) + {2-p)d{l + <^^(u))-P/^ 

and d > E{(p{u)); it implies X{u) <0ifp>2orp<2 and d < {p- l)/(2 - p). Henceforth 

O is increasing, and the same holds for P as a function of u* . Finally u* is decreasing with 
respect to /x, and consequently P is decreasing with respect to /z. □ 

Remark. When p = 2, and 17 = 2p — 1 = 3, equation (3.23 ) reduces to u" = —2u + 2u^, 
which, surprisingly, is an equation correponding to the problem with absorption, and (3.4 
) reduces to w" — w + = 0. In this case, all the solutions can be expressed in terms of 
elliptic integrals, see [3]. 
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4.4 Returning to the initial problem 

Proof of Theorem 2. Here [3 = (iq = p/ [q+l—p), A = A, is given by (1.15 ) and Cq = (3^~^Xq 
by (1.16 ). Moreover lu (a) = I3^q^w{j3qa) from (3.3 ), 6 = = -Cq/(3P and d = c/^p 



from (3.5 ). At end /(s) = g{s) = |s|' ^ s and h{s) = \s\'^ ^ s. Thus c> Cq is equivalent to 
b + d> 0, and then the constant sohitions w = ±{b + dY^^'^^P^^^ of (3.4 ) correspond to the 
constant solutions w = ± (c— Cq)^^^'^^^^^^ of equation (1.14 ). For any integer fc > 1, we 
look for periodic solutions lu of smallest period Itx jk. or equivalently solutions w of period 
Tfe = 2-K(3q/k. From (3.17 ), the function E is increasing. First consider the sign changing 
solutions: if c > Cq, then from Theorem 4.3, the period function T of w is decreasing from 
oo to 0, hence for any fc > 1 it takes precisely once the value Tk- If c < Cq, then T decreases 
from Td given by (4.4 ) to 0, thus it takes once the value Tk for any k > Mq = Td/lirfiq given 
at (1.18 ). Next consider the positive solutions: from Proposition 4.5, the period function 
of w takes any value between oo and l-nj ^{(i + 1 — p)(6 + d), thus it takes the value T^ for 
any fc < (p/3^~^(c — c,))^/^, which ends the proof. □ 

In the case of equation (1.14 ) (i.e. c = 0), we obtain the following description of the 
sets E and f +: 

Corollary 4.7 Assume p > 1, q > p — 1, and c = 0. 

(i) Then the set £ of changing sign solutions of (1.14 ) is given by (1.17 ), where kq ~ 1 if 
p<2 andq>2{p- l)/{2-p), and kg > M, if p > 2 or (p < 2 and q < 2{p - l)/{2-p)), 
where 

Mq = 2/{q-l), (4.10) 

if p = 2, and 



{{p - l)mq + l){m, - 1) ' ""''^ '"^ - V P{P - 1) ' ^ ^ ^ 



ifp + 2. 

If p > 2 or (p < 2 and q < 2{p - l)/(2 - p)), then £+ = 0. If p < 2 and q > 
2(p- l)/(2-p), then£+ = {(-Cg)i/(«+i-P)} . 



Proof. Here Cg < is equivalent to p < 2 and q > 2{p — l)/{2 — p). Furthermore Mq = 
To/27r/3g can be computed from (4.5 ), which gives (4.10 ), (4.11 ). Moreover in any case 
Cq+ Pq~^ /p = (iq{p " 1)(<Z + l)/p^ > thus there exist no positive nonconstant periodic 
solutions. □ 

Proof of Corollary 1. Let 5 be a sector on with opening angle 6 € (0, 27r). From [14, 
Th 3.3], ps is the positive solution of equation 

m) = + {^^s + - iPs + 1)') = 0, 

where k = n/O > 1. Using Corollary 4.7 (applied without assuming that fc is an integer) we 
distinguish two cases: 
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(i) p < 2 and q > 2{p — l)/{2 — p). Then there always exists a solution to the Dirichlet 
problem in S. Notice that < /3g < {2 — p)/{p — 1), thus ^(/Ss) < and consequently 

< Ps. 

(ii) p > 2 or p <2 and q < 2{p — l)/(2 — p). The existence is equivalent to A; > Mg (see 
(4.11 )). It means 

V fc; m,(p-2) y /32m2 ^^(^^ + _ 2)/3,/(p _ i)) ' 

Thus (?i(/3g) < 0. Equivalently, Pg < (is- □ 

5 The case p = \ 

5.1 Existence of a first integral 

As shown in Lemma 3.1, we can reduce the study to 

^ (^^=f ) - b + h{w) -d\w\-'w = 0, (5.12) 

where /i satisfies (3.8 ); in particular we are interessed by the case /i(s) = s. 

Here the problem is variational: if Si{w) is any primitive of w i— > \w^~^ fi{w) and 
R{w) if 6 7^ 0, i?(u;) = ln|u;| if 6 = 0, then (5.12 ) is the Euler equation of the 

functional 

n{w,w') = \w\''~'^ Vw'^ + - Si{w) + dR{w). 
Thus the following Painleve first integral is constant along the trajectories 

P{w, w') = /1\ - S^{w) + dR{w). (5.13) 

The system (3.9 ) reads as 
w' = y 

y' = G{w,y) = ^ — — > 

and it is singular on the line w = 0. For w > system (3.15 ) reduces to 

( W' = Wip{u) = Wy=f 

\ u' = bVT^^ - fi{w) + d. 
In the case f{w) = w, the equation satisfied by u is 



(5.14) 



u" = (1 - b)^=^u' -bu- d^=^. (5.15) 

V 1 — M V 1 — W 
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5.2 Existence of periodic solutions 



From the Painleve integral (5.13 ), we can describe the solutions, in the phase plane {w,y). 
Since a complete description is rather long, we reduce it to the research of periodic solutions. 



Proposition 5.1 Let p = 1, and consider equation (5.12 ). 

(i) Ifd^O, there is no periodic sign changing solution. Ifd = there exists such a solution 
if and only if b > —1, and then it is unique (up to a translation). 

(a) There exists periodic positive solutions if and only if b + d > 0. 

(Hi) Suppose moreover that fi{w) = w. Then the sign changing solution is given by 

w(t) = (6 + 1) cos(t — n); 

it has period 2tt. The orbits '7[(,j„o)] of the periodic solutions intersect the axis y = at a 
first point (/i, 0) such that ^ < a = b + d, and fi describes jjL € a) with p = ifd<0, and 
fl> if d> 0; it is given by (5.20 ),(5.18 ),(5.19 ). 

Proof. By symmetry we reduce the study to the case > and the painleve integral (5.13 
) takes the form 

Vl - - Si{w) + dR{w) = C, (5.16) 



w 



where we denote 



/ s''-^f^{s)ds 
Jo 



and 



Si{w)= s''-'^ fi{s)ds if6>-l, 



s*-Vi(s)rfs+^ if6<-l, 

/W 
s-^fi{s)ds if6 = -l. 



Step 1. Periodic sign changing solutions. The curves in the phase plane {w,y) are given, 
for w > 0, by 



2 _ / ^'^^ \ .,,2 



^ " \C - dR{w) + Si{w) ' ~^ 

~ {Si{w) + C - dR{w)f ' 

which defines zby in function of w. If there exists a sign changing periodic solution, the 
trajectory intersects the axis w = at some point (0,£) with ^ > 0, thus y needs to ^ as 
w tends to 0. Prom (5.16 ), it is impossible if 6 < —1. Assume 7^ 0; if — 1 < b, then near 
w = 0, in any case < {b"^ /d^ + l)w^, thus 1 = and w' /w is bounded, thus the maximal 
interval of existence is infinite, and we reach a contradiction. If d = 0, and C ^ 0, then 
y"^ = -w'^{l + 0(1)), which is impossible. If = C = 0, then = w'^{w'^'' /Sf{w) - 1) ; 
observing that the fimction w 1— > x(^) ^ w~'^Si{w) is increasing from to 00, the curve 
intersects the two axis at (0, 6+1) and (x~^(l), 0) and this corresponds to a closed orbit. 
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Step 2. Existence of periodic positive solutions. If we look at the intersection points of any 
trajectory in the phase plane with the axis y = 0, we find that they are given by H{w) — C, 
where 

H{w) =w'' + dR{w) - Si{w). 

Then H'{w) = w^~^{b + d — fi{w)). If 6 + d < 0, then H is decreasing, thus there ex- 
ist no positive periodic solutions. U b + d > 0, the function H is increasing on (0, a) 
where it reaches a maximum M, and decreasing on (a, oo). The stationary point (a, 0) with 
a = fi^{b + d) corresponds to C = M. If 6 > 0, then lim^u^o-ff = 0, while, if 6 < 0, 
then lim^^o-ff = — oo. Equation H{w) = C has two roots < «;i < W2, if and only if 
C e (maxjlim^^o H, lim^^ao H} , M). Moreover, if there exists a trajectory going through 
(wi,0) and (w2,0) and if one denotes 

K{w) = dR{w) - Si{w) = H{w) - w'', 

one has K{w) < C on {wi,W2), thus C > M' = max K = K{f^^{d)). Conversely, if 

maxjlim ff, lim ff,M'| <C<M, M = H{f-\b + d)), M' = K{f-\d)), (5.17) 

then there exists a closed orbit going through {wi,0) and {w2,0). 

Step 3. End of the proof. The sign changing solution is given by u>^ + w'^ = (6 + 1)^ and 
its trajectory is a circle with center and radius b + 1; for w > 0, w = {b + 1)^1 — = 
by/1 — — u', thus u' = — \/l — u'^, and 6' = —1, then w{t) = (6 + 1) cos(^ — ri), periodic 
solution with period 27r. Now consider the positive periodic solutions. Here a = b + d, and 

' {l + d/b)w^ -w^+'^/{b + l), if 6 7^ 0,-1, 
H{w)=i 1 + dlnw-w, if 6 = 0, (5.18) 

(l-d)w-i -Inw, if6 = -l, 

M = a^+^/b{b+l), M' = {d+f+^/b{b+l)), if 6 ^0,-1, 

M = l + dhid- d, M' = dlnd-d, if 6 = 0. (5.19) 

M = -l-ln{d-l), M' = -1-lnd)), if6=-l. 

If d < 0, thus b > 0, then any C € (0,M) corresponds to a closed orbit, thus for any 
fi & (0, ai) , one has a closed orbit passing through {n, 0) , of period still denoted by T+(^). 
If d > 0, in any case, any C € (M',M) corresponds to a closed orbit. If — 1 < 6 < 0, 
then H is increasing on (0, a) from — oo to M < 0, and then decreasing on (a, oo) from 
M to -00. If & < -1, then M > M' > 0. If b < -1, then d > 1, and limt„^o ^ = -oo, 
lim^_>oo -ff = 0, < M' < M. If 6 = 0, thus d > 0, then hm^^.^o H = — oo, then any 
C e (M — 1, M) corresponds to a closed orbit. Then H is increasing on (0, d) from — oo to 
M = 1 + dhid — d > 0, (notice that M — d = 1) and then decreasing on (ai, oo) from 
M to — oo; let p. G {Q,b + d) be defined by 

H{fl) = M', (5.20) 
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thus for any ^ e {fl, a) , one has a closed orbit passing through (/x, 0) , with a period 
still denoted hy T+{fi). If -1 < b < 0, thus d > -b > 0, then lim^u^o = -oo = 
lim^^oo H, M < 0, and any C E {M' , M) corresponds to a closed orbit (if b = —1, then 
H{w) = (1 — d)w~^ — \nw, M = —1 — ln{d — 1), M' = — 1 — \nd). Returning to equation 
(1.14 ), the conclusion follows with jig = jl''. □ 



5.3 Period of the solutions 

Let p=l,b + d>0. Consider the equation (5.12 ). Let T'^{fi) be the least period of the 
periodic positive solutions corresponding to the orbit 7[(^,o)]- As in the case p > 1, we have 
a general result: 

limr+(M) = ^_. (5.21) 

Next we study the variations of the period in the case of a power fi{w) = tu. 

Theorem 5.2 Assume p = \,b + d> Q and fi{w) = w. Then lini^^a T^{ji) = lixj^Jb + d. 
If d < 0, then lim^^o T'"'" (m) = oo. If d > 0, then lim^^pT+(^) — r+ is finite, and given 
by (5.23 )ifbi {0,-1}, by (5.24 ) if b = 0, and by (5.25 ) if b = -1. If d = 0, then 
T+ = tt{1 + 1/6). 

Proof. Step 1. Assume b ^ {0, —1}. From (5.16 ), the solutions of (5.12 ) satisfy 

d 



6 + 16 



thus 

= 6^1^ - «; + d = -^/^^ - ^ + 

6 w" 

Eliminating w between the two relations, we find that C6(6 + 1) > and 

(d + bVl - V? - u') '^^'^'^ (d + 6^1 - V? + bu'') '^^''^'^ = (C6(l + 6))i/(^+i) := A. 

When a solution goes through the half-part of its trajectory T located in Q, the associated 
function u increases from to some u* € (0, 1) where the derivative u' vanishes and d + 
6Vl - > 0; next d+bVT — V? is monotone and positive at and u* , thus d+b^/1 — > 
everywhere. And A = d + b^/l — w*^ = w* (the value of w when u = u*). Let 

and G{s) = {l-sf^''+'Hi + bsf^''^'K 



d+bVT^ 



If 6 > 0, then.2 e (-1/6, 1); if 6 < -1 then2; e (-oo, 1) ;if-l < 6 < Othen.2 e (-oo, 1/|6|) , 
and 

G{z) = ^ 

^ ' d+bVT^ 

Since 

G'{s) = -bs (1 - 8)-'^^'+'^ (1 + 6s)-^/(''+^) , 
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and 

G"{s) = -b{l- (1 + 5,)-(2f+i)/(6+i) ^ 

it follows G(0) = 1, and is a maximum if 6 > and a minimum if 6 < 0: if 6 > 0, G 
increases on (— 1/&, 0) from to 1 and decreases on (0, 1) from 1 to 0; if 6 < 0, G decreases 
on (—00,0) from oo to 1 and increases on (0,min(l, 1/ |&|)) from 1 to oo. Thus it has two 
inverse functions —Li and L2 : for 6 > 0, Li maps (0, 1) into (0, 1/6) and L2 maps (0, 1) 
into (0, 1) ; for 6 < 0, Li maps (1, 00) into (0, 00) and L2 maps (1, 00) into (0, min(l, 1/ |6|)). 
Then ^ 

T+ = r+ + T+, T+= [ V-i.u. (A)dA, (5.22) 

Jo 

where 

'Wz,u' (A) = 



{d + bVl - X^u*'^)Li{{d + bVl-u*'^)/{d + bVl - X^u*'^)) ' 

• First suppose d < (thus 6 > 0); then one looks at the case where G — > 0, thus 
Vl - w*^ -d/b, thus u* -*u= ^Jl^^'WJW- . Near u, 



b{Vl - A2u*2 - VI - M*2)L,(0) - 6Li(0)(l - A2) ' 

therefore tends to 00. 

• Suppose d > 0,b > 0. One looks at the case where G ^ M', thus m* ^ 1. There 
exists a constant ni > such that < 1 — G(.s) = G(0) — G(s) < m^s^ on [—1/6,1]. 
Indeed G'(0) = and G" is bomided on [-1/26.1/2], and on [-1/6,-1/26] U [1/2.1] the 
quotient (G(0) — G{s))/s^ is bounded. Thus 1/Li{ri) < m/y^l — rj on [0, 1) , hence taking 
rj = (d + 6\/l — u*'^)/{d + by/1 — A^u*^), and computing + 

6(1 - A2)u*2 

1 — T] = 



{d + bVl - A2m*2) (VI - U*2 + VI - A2u*2) 



one finds Vi,M* (A) < Amj yjb (1 — A2). Prom the Lebesgue theorem, as u* ^ 1, T+ tends to 
the finite limit 

y+ _ I 2^+ 2=-+ = 2 /" ^ ^^^^ (5 23) 

' ^' ' io (d+6vr^)ii(rf/(d+6\/r^)) 

in particular if = 0, then Li(0) = 1/6, Li(0) = 1, thus T+i = tt and f+j = 7r/6. 

• Suppose 6 < 0, thus d > — 6 > 0. Then again G M', consequently u* ^ 1. The 
function 

u* ^Q{u*,X) = r]= - — -^==== = 1- 



d + 6V1 - A2u*2 ((^ + 6V1 - A2u*2) (Vl - u*2 + Vl - A2m*2) 



is increasing on (0, 1) from 1 to d/(d+6Vl - A2) and d/(d+6Vl - A2) < d/{d+b) = a. There 
exists m > such that < G(s) — 1 < m2s2 on [— Li(a), L2(a)] , thus 1/Li{r]) < m/\/r] — 1 
on {l/d/{d + 6)] . Thus as above, ipi^u* (A) < ^m/-s/\b\ (1 - A2), and r+ tends to 7"+ defined 
at (5.23 ). 
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Step 2. Assume 6 = 0. There exist periodic solutions for any C G {M — 1, M) . The solutions 
are given by 

\/l — u'^ + H{w) = — V? + d\nw — w = C 
and u' = —w + d, thus u is maximal (= u*) for w = d : therefore \/l — u*^ + H{d) = C, then 



H{d - u') = H{d) + Vl - w*2 - Vl -u2 

and has two inverse functions Hi from {—oo,H{d)) into (0,d) and (cJ, oo), thus (5.22 ) 
holds with 

' (d - H,{H{d) + VI - u*'^ - Vl - A2u*2) 



and ^ = + Vl - - Vl - A2m*2 = if (d) - A: = + h) stays in (M - 1, M) = 

(il(rf) - 1, H{d)) , and H{d+h)-H{d) > -m^h'^ for e (M - 1, M) , thus H{d)-^ = 

k < m?{d-Hi{^)f, thus 

2m 2m (Vl - w*^ + Vl - A2u*2) 4m 
V'i.u* (A) < ^ = ^ , ^ < 



Vfc vn^ " \/r^rx2' 

Therefore, as m* ^ 1, T+ tends to the finite limit 

d\ 

{d-HiiH'idj^Wf^)' 



f+ = f+ + f+, i;+ = 2^ — 



(5.24) 



Step 3. Assume b = -1. In that case d > 1; let B = -(C + 1) G (ln(ii- l),lnc!) then 
B — > In and 

u' + w = d — \/\ — V? = (B + \)w — wlnw = Hb{w) 

where Hb is increasing on (O, e^) from to and decreasing on (e^,oo) from e^ to 
— oo; it has two inverse functions LB,i from (— oo,e^) into (0, e^) and (e^~^,oo); and 
w* =d- Vl - u*^ = e^; then (5.22 ) holds with 

- VI - A2«*2 - LB,i{d - VI - A^M*^ |iis-i(is,i(rf- VI- A^u*^))!' 

Because Hs-iie^) = 0, Hb~i{x) — HB-i{e^) — H'g_-^{(,){x — e^) and x ranges onto 
{HB,i{d- l),HB,2{d- 1)) (a:i,B, a;2,B), when B Ind, (xi^s, 2:2,3) ^ (a:i,ind, a;2,in<i), 
it follows \H'g_^{^)\ > l/jjL > independent on B. Moreover iJB(a;)-ilB(e^) = {l/2)H'^{S,){x- 
e^)2 = -(l/2^)(x - e^)2. Thus there exists m > such that 

Hb{x) - HB{e^) < m^{x - e^f < HI_^{x). 

Therefore, near Ind, taking x = LB,i{d — Vl — X^u*^), one derive 

2 2 4 

V" ♦ W < 1=^^^^^^^= = : < . 

m/uV - Vl - A2u*2 - m/xV Vl - w*^ - Vl - X^u*"^ m/iVl - A2 
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Consequently, as u* 1, tends to the finite limit 



T+ = T+ + r+ n, = ^fT7; —^^—^==-. (5.25) 

Jo |-Hlnd-l(i'lnd,i(ct- Vl - A^))| 

□ 

Remark. In the case = 0, 6 ^ 1, notice that and T2 converges to 7r/\/6 as jj, tends 
to h (one can verify it by linearizing the equation in u) and respectively to tt and ir/b as ^ 
tends to Thus if those functions are monotonous, they vary in opposite senses and it is not 
easy to get the sense of variations of their sum T+. Moreover in the phase plane {w,y\), 
as /i tends to 0, on can observe that the trajectory tends to a limit curve constituted of a 
segment [(0, 0), (0, 6)] and half of the unique closed orbit surrounding (0, 0), circle of center 
and radius b + 1, which is covered in a time tt 

The case 6 = 1 is the most interesting for (5.12 ), since it corresponds to the initial 
problem (1.14 ). In that case we improve the results by showing the monotonicity of the 
period function: 

Theorem 5.3 Assume b = 1, d > —1. When d = the period function T'+(/u) is constant, 
with value 2tt, thus there exists an infinity of positive solutions w of (5.12 ), which are all 
2n-peHodic; they are explicitely given by 



w 



= Vl-K*'^ sin^ T-K* cos r, re[-7r,7r], K*g{0,1). (5.26) 



When d ^ 0, then T~^{fi) is strictly monotone; if d < it decreases from 00 to 2^/^/1 + d; 
if d > it increases from 

^{d+Vl^f^^ V COS0 + 2d 

to 27r/\/l + d. 

Proof. • If d = 0, then u" = —u, from (5.15 ), and u = sin^ G [0, 1) , thus the positive 
solutions w are given in Q by 

u = K*smT, K*€[0,1), Te[0,7r], 

and the period T+ is constant, equal to 27r. We obtain an infinity of positive solutions w, 
given explicitely by 



w 



= \J\-u^- u' = Vl - sin^ T-K* cos r, K* G (0,1) 



which intersect the axe y = at points = (1 =F K*). 
• In the general case d> —1, we find 



[d + V1-m2 - u'^ (v'1-m2 +u' + d}j=A^ 
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that means G is symmetric: G(s) = \/l — s^, thus 

u'^ = {d+^/l^f-A^ 
a/I - u*2 = A-d = \/2C - d; thus here = , and 



^+ , f dX ^ . ^. {d+s/l^^sf-{d+^/T^sf 
=4 , -., where *(s, A) = -^^ ^ —. 

We show that the period function is strictly monotone with respect to u* . Because 

s2^(s, A) = d{d + 1) (l/Vl^s - 1/Vl - A2s) >0, 

we see that T+ is increasing if d < and decreasing if d > (and we find again that it is 
constant if = 0). Also fi can be expressed explicitely in terms of u* by 



^ = d+l-\l{d+lY-{d+\/l- M*2) , 



Therefore /i is decreasing with respect to u* , hence T+ is decreasing with respect to if 
d < and increasing \i d> Q. □ 



5.4 Returning to the initial problem 

Proof of Theorem 3. Here (iq = I3q = 1/q, the substitution (3.6 ) takes the form uj {a) = 
\w{a)\^^'^ ^ w{(7), and thus 6=1, and d = c from (3.7 ). Then the existence of sign changing 
solutions of (1.14 ) is given by Proposition 5.1. The constant solutions exist whenever 
c+ 1 > 0. Next we look for positive solutions of smallest period 27r/A; applying Theorem 5.2 
and Theorem 5.3. If c < the period function T+ decreases from oo to 27r/-\/l + c > 2tt, 
thus there exist no solution. If c > 0, r+ increases from r+ given by (5.27 ) to 2^^/^/T+^, 
thus it takes once any intermediate value, which gives one solution (up to a translation) 
for any k G {ki,k2)- If c = 0, the solutions uik are given explicitely by (5.26 ), and u)q is 
obtained from loq; this means that system (3.9 ) does not satisfy the uniqueness property at 
(0,0). □ 
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